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Abstract. We present a new tool for the photometric estimate of stellar masses in distant galaxies. The observed 
source spectral energy distributions are fitted by combining sets of various single stellar populations, with different 
normalizations and different amounts of dust extinction, for a given (Salpeter) IMF. This treatment gives us the 
best flexibility and robustness when dealing with the widest variety of physical situation for the target galaxies, 
including inactive spheroidal and active starburst systems. We tested the code on three classes of sources: complete 
samples of dusty ISO-selected starbursts and of K-band selected ellipticals and SOs in the HDF South, and a 
representative sample of z ~ 2 to 3 Lyman-break galaxies in the HDF North. We pay particular attention in 
evaluating the uncertainties in the stellar mass estimate, due to degeneracies in the physical parameters, different 
star formation histories and different metallicities. Based on optical-NIR photometric data, the stellar masses are 
found to have overall uncertainties of a factor of ~ 2 for E/SOs, while for the starburst population these rise to 
^— ^ ' factors 2—5 (even including ISO/15/im photometric data), and up to > 10 for Ly-break galaxies. Our analysis 

' reveals in any case the latter to correspond to a galaxy population significantly less massive (M < a few 10^°Mq) 

, than those observed at lower redshifts (for which typically M > several 10^" M©), possibly indicating substantial 

stellar build-up happening at 2 ~ 1 to 2 in the field galaxy population. Using simulated deep SIRTF/IRAC 
observations of starbursts and Lyman-break galaxies, we investigate how an extension of the wavelength dynamic 
range will decrease the uncertainties in the stellar mass estimate, and find that they will reduce for both classes 
O ' to factors of 2—3, comparable to what found for E/SQs and good enough for statistically reliable determinations 

^ , of the galaxy evolutionary mass functions. 
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1. Introduction the true rates of star formation are affected by a variety of 

uncertainties and biases, expecially related to the amount 
One of the still open critical questions of modern cos- ^^^^ galaxies and its effect on the SFR tracers, 
mology is to understand the epoch at which galaxies as- 
sembled the bulk of their stellar content. In the so-called A complementary approach consists in measuring the 
monolithic scenario, the assembly of galaxies took place dynamical or stellar masses of distant galaxies, instead 
on rapid timescales at high redshifts, then galaxies evolved of the instantaneous SFR. Dynamical masses are directly 
passively to present days. On the contrary, in the hierar- connected to theoretical predictions, but very difficult 
chical scenario (Kaufmann & Chariot, 1998) galaxy for- to measure, requiring high spatial and spectral resolu- 
mation is predicted to be a more continuous process and tion spectroscopy of selected samples of faint high-redshift 
elliptical galaxies to assemble through merging of lower galaxies. Stellar masses, on the contrary, are less well de- 
mass disc galaxies at moderate redshifts. termined by theory, but have the advantage that can be 
Many authors have analysed the evolution of the global derived using multiband optical and near-IR photome- 
star formation rate (SFR) through cosmic history, back try as a powerful alternative to time-expensive spectro- 
to z>3, using different observational tracers of the SFR scopic investigations. Indeed galaxy near-IR SEDs show a 
in distant galaxies. It is generally accepted that galaxies moderate dependence on the age of the contributing stars 
produced stars more actively in the past than today, but (e.g. Franceschini & Lonsdale, 2003), because in a typical 
galaxy the stellar mass is dominated by low-mass stars, 
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1999, Origlia & Oliva 2000), these moderate-mass stars 
emit predominantly in the near infrared (NIR, J-to-K rest- 
frame bands) , especially during their cool giant phase, and 
are only slightly affected by dust extinction. 

In this work we present a new spectro-photometric 
synthesis tool, aimed at the estimate of stellar masses in 
galaxies from a thorough analysis of their spectral energy 
distributions (SEDs). The tool combines a set of single 
stellar populations of different ages, assuming different 
star formation rates and different amounts of extinction 
for each. The code is tested by fitting the optical-NIR 
broadband SEDs of two samples of intermediate (z ~ 0.5 — 
1.5) redshift galaxies, luminous IR-selected starbursts in 
the HDFS (Pranceschini et al., 2003) and mostly passively- 
evolving K-band selected ellipticals from Rodighiero et al. 
(2001). A benchmark higher-redshift population, also tar- 
geted by our analysis, was extracted from a sample of 
2: ~ 2 — 3 Lyman-break galaxies (Papovich et al., 2001). 

We pay particular attention to the uncertainties in the 
mass estimates, due to degeneracies in extinction, age and 
star formation history, and discuss how forthcoming near- 
to-mid infrared data from SIRTF will further constrain 
the photometric mass estimate for these sources. 

The paper is structured as follows. Section 2 describes 
the three samples of selected sources. Section 3 presents 
the model tools used to fit the observed optical-NIR spec- 
tral energy distributions of the sources and the different 
star formation histories assumed for the two classes of 
galaxies analyzed. Section 4 discusses our results, with a 
particular care for degeneracies. Section 5 shows perspec- 
tives for future SIRTF/IRAC observations and Section 6 
summarizes our conclusions. We assume a ifg = 65 [km 
s~^ Mpc~^], = 0.3, Oa = 0.7 cosmology. 



2. The samples 

Our analysis compares synthetic spectra of a new spec- 
trophotometric code with broad-band observational data 
to estimate the galaxy's main physical parameters, partic- 
ularly the total stellar mass, but also, at the same time, 
age, extinction and star-formation rate. Our aim here is 
to investigate in particular the uncertainties in the stellar 
mass estimates from degeneracies in the solutions due to 
the sparse characterization of the observational spectrum. 

We compare with test cases provided by populations 
of luminous starburst galaxies and spheroidal galaxies, the 
two representing somewhat limiting physical situations. 
In the former class, the presence of red giants and dust- 
reddened young stars complicates the analysis and widens 
the uncertainty range. On the other hand, for passively 
evolving spheroids the uncertainties are expected to be 
much lower. 

A third situation explored here is higher redshift {z ~ 
2 — 3) Lyman-break galaxies, to show how the poorer cov- 
erage of the rest-frame spectral energy distribution is af- 
fecting the uncertainties in the parameters. 



2.1. Mid-IR selected starbursts 

The rich variety of high photometric quality data in the 
Hubble Deep Field South and its Flanking Fields allows 
us to build accurate SEDs for galaxies. We have selected a 
sample of intermediate-redshift sources detected by ISO 
in the LW3 (Ac = 15 /xm) band (Franceschini et al. 2003, 
Oliver et al. 2002) to a flux Hmit of Slws = 90 /iJy. 
Franceschini et al. (2003) compare the ISO source list with 
those from the Deep ESO Imaging Survey (EIS Deep, da 
Costa ct al. 1998) and with the optical UBVRI catalogue 
by Teplitz et al. (1998), and identify 35 extragalactic ob- 
jects. Of these 35 sources, 3 are likely to be hosting an 
AGN component and for another 2 not enough photo- 
metric data are available. For the remaining 30 galaxies, 
UBVRIJHK photometry is available (Da Costa et al. 1998, 
Tepliz et al. 1998), as well as the LW3 15 /xm ISO flux and 
— in some cases — the LW2 6.7 /um flux. For 17 sources 
spectroscopic redshifts are available from Rigopoulou et 
al. (2000) and Franceschini et al. (2003); for the 13 others 
we rely on the photometric estimate. 

2.2. K-band selected ellipticals and SOs 

As representative of this class we chose a morphologically 
selected E/SO sample by Rodighiero et al. (2001) in the 
HDFS. Based on a SExtractor analysis (Bertin & Arnouts, 
1996) of the EIS-Deep K band SOFI image (Da Costa et 
al., 1998), they computed a completeness limit of K~20.2 
(Vega magnitude). Among this list, the subpopulation of 
spheroidal galaxies was identified from both visual mor- 
phological inspection and using the authomatic analysis 
tool GASPHOT (Pignatelh & Fasano, 1999). From this, a 
sample of 29 E/SO brighter than K=20.1, in the WFPC2 
HDFS field was selected. 

Among these sources, data in the whole UV-to- 
NIR spectral range are only available for 23 galax- 
ies: UBVI from HST (filters F300W, F450W, F606W, 
F814W, Williams et al. 1996) and JHK from EIS Deep 
(SOFI/NTT, Da Costa et al. 1998) photometry. For 10 
of these sources spectroscopic redshifts arc available from 
Vanzella et al. (2003) and Sawicki & Mallen-Ornelas 
(2003). For the remaining, the photometric estimate is re- 
liable to within 5%. 

2.3. Lyman break galaxies 

Papovich et al. (2001) report the analysis of stellar popula- 
tions of a sample of Lyman-break galaxies in the Hubble 
Deep Field North, selected through the U dropout tec- 
nique, up to redshift 3.4. 

The publicly available photometric data consist of 
the original optical WFPC2 F303W, F450W, F606W, 
F814W HDFN observations by Williams et al. (1996), 
the NICMOS F160W and ground based (IRIM) Ks data 
by Papovich et al. (2001) and the J band magnitudes re- 
ported by Sawicki & Yee (1998), obtained with IRIM on 
the KPNO 4m Mayall telescope. Redshifts have been spec- 
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Fig. 1. HDFS/LW3 source S40 ai z = 1.27: best fit mass is M = 1.63ell. Left Panel: comparison of the observed 
SEDs to the best fits. The sohd hne is the sum of all the used SSPs; the dashed hue is the sum of aU the young SSPs 
(age < 10^ yrs); the dotted line represents the old populations (age > 10^ yrs); the dashed thick line longward 5 /Ltm 
(restframe) is the k-corrected rescaled M82 template; the square black dots are the observed photometric datapoints; 
the triangular red point at A ~ 15 is the synthetic LW3 flux (overlapping to the observed LW3 when the FIR 
constraint is satisfied). Right Panel top: values of SFR (histogram, left y-scale) and E(B-V) (dots, right y-scale) for 
each single stellar population. Middle: ratio between the FIR and V-band fluxes for the different SSPs. This quantity is 
used to asses both the contribution of various stellar populations to the FIR luminosity and the effect of extinction on 
each population (depressing the optical light and boosting the IR). The global ratio between the FIR and V fluxes of 
the composite final synthetic spectrum is reported on the right. Bottom: cumulative stellar mass formed in the galaxy. 
The youngest stars, significantly contributing to the FIR bolometric emission, provide almost no contribution to the 
total luminous mass assembled in the galaxy. 



troscopically derived by various authors (e.g. Steidel et al. 
1996, Lowental et al. 1997). 

We have selected a sub-sample of 8 sources between 
z ^ 2 — 5: this redshift range is the one where future 
SIRTF/IRAC observations will be more effective in con- 
straining the stellar masses of such objects. 



3. The Model Tool 

The observed SEDs for the target galaxies span wide 
ranges of rest-frame wavelengths (0.2—1.4 /im at the mean 
redshift of the starburst sample, z = 0.6), and therefore 
are sensitive to stellar populations over a wide range of 
ages and mass/luminosity ratios and to different aspects 
of galaxies star formation histories. In the starburst galaxy 
case the observed UV-optical SED is strongly affected by 
dust extinction. Dust re-radiation is expected in the mid- 
to far-IR. The relative contributions of young red stars 
from older stellar populations can then be constrained 
with ISO mid-IR observations. 



3.1. Generalities on the spectro-photometric model 

To generate the synthetic spectra, we have developed a 
new automatic analysis tool as an upgrade of the spectral 
synthesis code by Poggianti et al. (2001). We modified 
it on the computational side, in versatility and improved 
user-interface. 

Each observed SED is modelled as a combination of 
a set of single stellar populations (SSP) of solar metallic- 
ity and different age. Each SSP is meant to represent a 
formation episode of average constant star formation rate 
(SFR) over a suitable time period At. Each single popu- 
lation is assumed to be extinguished by dust in a uniform 
screen according to the standard extinction law of the dif- 
fuse medium in our Galaxy {Ry = Av/E{B — V) — 3.1, 
Cardelli et al. 1989; see Poggianti et al. 2001 and Berta 
et al. 2003 for a discussion). The extinction value E(B-V) 
is allowed to vary from one stellar population to another, 
and the total spectrum is built up by summing the extin- 
guished spectral energy distributions of all generations. 
However, in consideration of the fact that high extinction 
values are found only for the young stellar populations 
embedded in thick molecular clouds and that disc popula- 
tions are on average affected by a moderate Ay (< 1, e.g. 
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Fig. 2. Left: contours on the M vs. Lpjn diagram, for source S40. Contours are drawn at 1,2,3a levels, corresponding 
to 68.3%, 95.5% and 99.7% confidence (decreasing as shade darkness increases). The solid and dashed horizontal lines 
represent the far-IR huninosity and uncertainty as derived from the observed LW3 flux, based on fitting with the M82 
template. Right: values as a function of the stellar mass of the models, for those solutions reproducing the LW3 
flux within 1%. The mass uncertanty range is [0.70,3.80] lO^^M©. 



Kcnnicutt, 1992), wc have limited the E(B V) values for 
the populations with ages > 10^ yrs to be less than 0.30. 

All the SSP spectra have been computed with a 
Salpcter initial mass function (IMF) between 0.15 and 120 
Mq, adopting the Pickles (1998) spectral library extended 
with Kurucz (1993) atmosphere models. The composite 
spectra, including both photosphcric stellar and nebu- 
lar (line and continuum) emission, have been obtained 
through the ionization code CLOUDY (Ferland, 1990). 
This procedure provides a reliable description of single 
stellar generations up to a wavelength of ~ 5 /xm (rest- 
frame). Beyond this wavelength, dust emission becomes 
typically no longer negligible. 

The bulk of our analysis assumes SSPs of solar metal- 
licity. However we can also deal with synthetic spectra 
with different metallicities, i.e. one third and one tenth 
solar, Z=0.008 and Z=0.002. As it wiU be subsequently 
discussed, the estimate of the stellar mass content of a 
high-redshift galaxy barely depends on the detailed as- 
sumption for metallicity. 

The model SED is obtained by combining a number 
Nssp of different SSPs. weighted by different values of 
SFR and different amounts of extinction: 



j^unext ^ j^SSP . gpj^, . 



mod 



Nssp 

E 



10 



(1) 



(2) 



The sum is intended over all the Nssp single stellar gen- 

3; 



erations; L'^".'^^* and Av,i are respectively the unextincted 



spectrum weighted by SFRi and the extinction value for 
the i-th SSP, if f ^ is the unextinguished luminosity of 
the i-th SSP normalized to 1 M©, Ati is the duration of 
the i-th SSP and finally ^ is the standard extinction law 
{Rv = 3.1). 

The result is corrected for cosmological dimming by 
applying the K-correction and compared to the observed 

SED with: 



K{v,z) 



K{u,z) = {l+zy''' ^ 



z)]T{v)dv 



CL{y)T{u)dv 



(3) 



(4) 



where is the luminosity distance at the given redshift 
(Carroll et al., 1992) and K{v, z) the K-correction factor; 
T{v) is the filter transmission function. 

3.2. Fitting the starburst galaxy SEDs 

The SEDs of the sources in our mid-IR sample have 
been modelled as a combination of up to 10 SSPs, whose 
ages have been chosen by considering the evolutionary 
timescales of the associated stars. We have then con- 
sidered up to four young solar-metallicity stellar gener- 
ations (with ages of lO'', 3 lO'', 8 lO'', lO'^ yr), mostly re- 
sponsible for the ionizing photons producing the emission 
lines and for dust heating and FIR emission; up to five 
intermediate-age populations (5 lO'^, 10^, 3 10^, 5 10^, 10^ 
yr) with the strongest Balmer absorption lines; finally for 
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Fig. 3. Comparison of three different fits of the SED of source E5 in the elhpticals HDFS sample by Rodighiero et al. 
(2001), obtained with the three different adopted methods. Left: Analytic star formation history, given by a Schmidt 
law; center: two single stellar populations (10^ and 1.2 10^° yrs old) with solar metallicity; right: two old populations 
(1.2 IQi" yrs) with different metallicities (Z=0.02, 0.008). TOP: the best fit models SEDs (solid hne) have been spht 
into two components: on the left into the contributions of SSPs younger than 10^ yrs (dashed line) and older (dotted 
line); at center into the contribution of the 10^ yrs old SSP (dashed) and of the 1.2 10^" yrs old one (dotted): on the 
right into the contribution of the Z=0.008 SSP (dashed) and of the Z=0.02 one (dotted). BOTTOM: star formation 
history and contribution of each population to the total stellar mass assembled in the galaxy. In the bottom-right 
panel the two old SSPs with different metallicities have been distinguished both in the SFH plot (shaded histogram is 
the Z=0.02 SSP) and cumulative-mass diagram. The main contributions to the for the central and right solutions 
come from the J and V data respectively. No extinction is plotted, since it was assumed constant for all populations 
(see Table 



older generations of stars we have included a population 
modelled with a constant star formation rate (SFR) be- 
tween 2 and 12 Gyr before the observation. This latter 
choice is discussed in section l?!^ together with the effects 
of varying the SSP's metal content. 

In the attempt of reducing the uncertainties in the 
mass estimate due to the contribution of young stars to 
the near-IR restframe SED, we have considered the ISO 
15 /zm flux in addition to the UBVRIJHK photometric 
data. The dust emission spectrum longwards of 5 fiia 
(restframe) is modelled by assuming the observed M82 
IR SED and then comparing to the flux detected by ISO 
in the LW3 band (12 to 18 ^m). The M82 IR template is 
adopted as a typical spectrum of luminous IR starbursts 
(except some exceptional ultraluminous galaxies like Arp 
220). The predicted IR emission is calculated as the dif- 



ference between the unextinguished and the extinguished 
optical spectrum, assuming that all the flux absorbed by 
dust is reprocessed and re-emitted longwards of 8 /xm. 
The M82 template is then rescaled in such a way that 
its bolometric emission between 8 and 1000 /im coincides 
with the dust reprocessed luminosity. Finally, the prop- 
erly scaled M82 SED is K-corrected and convolved with 
the LW3 filter reponse to compare with the observed LW3 
flux. 

The M82 template was adopted to model the IR SED 
of the ISOCAM sources as a fairly conservative choice. 
Elbaz et al. (2002) have shown that tight nearly linear cor- 
relations between the observed mid-IR emission and the 
bolometric FIR luminosity of local IR galaxies hold over 
^ 4 orders of magnitude in Lju (between 10^ and 10"'^^ 
Lq), involving IRAS 12 /Ltm and ISO LW2, LW3 observed 
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Fig. 4. Best-fits spectral solutions for representative elliptical galaxies, based on the Schimdt star formation history 
model (see Sect. l3.3fl . See also Figs.n]and|21 



fluxes. M82, with Lm ~ 3.3 10^° Lq, is a prototypical 
galaxy, representative of dusty starbursts and laying in 
the middle of the Elbaz et al. (2002) correlations. 

Objects not following these relations, e.g. the ULIRG 
Arp 220, are characterized by a much more peaked FIR 
emission and dust self-absorption. Adopting such an ex- 
treme galaxy as IR template would produce higher bolo- 
metric infrared luminosities for a given 15 /im flux, and 
therefore require very high — often unphysical — extinc- 
tion values for those SSPs contributing to the FIR emis- 
sion. The estimate of stellar mass and its uncertainty (see 
Sect. I4.2|l is in any case not significantly affected by this 
choice. 

The exploration of the parameter space (see more de- 
tails in Sect. 3.5 below) starts with a large number (10) of 
stellar populations. However, after the minimization pro- 
cess and for each solution with acceptable x^, SSPs provid- 
ing less than 1% of the observed spectrum at each wave- 
lengths (after extinction has been applied) are identified 
and excluded from further consideration. 

The effective number of SSPs involved in the best-fit 
solution of each examined object is reported in Tabled 
Tipically 3—4 stellar populations then contribute to the 
fit of the observed SEDs, which in turn means a total of 
6—8 parameters (SFR and E[B-V] for each population). 

The optimization process involves, at least at the be- 
ginning of the exploration, a large number of popula- 



tions/parameters to ensure that the most general descrip- 
tion of galaxy spectra and the most conservative evalua- 
tion is given of the uncertainties in the physical parame- 
ters. 

3.3. Fitting the elliptical SEDs 

It is common wisdom that the properties of local ellipti- 
cals/SOs may be reproduced by one suitable Simple Stellar 
Population with old age. However in some cases there are 
hints for the presence oi a) a, younger population, repre- 
sentative of a partial renjuvenation episode occured in the 
recent past; h) oi a significant contribution from sub-solar 
metallicity stars. Both effects result in a bluening of colors. 

We have analysed the spheroidal galaxy population 
with various methods. The first approach that we have 
considered consisted in adopting a simple analytical form 
of the star formation history (SFH) in E-type galaxies. 
We have achieved this by combining 55 SSPs in an almost 
continuous age sequence in which the stellar mass formed 
in each generation is obtained from: 

s™,«.(ii^)'%cxp(-it_i) (5, 

where is the epoch of formation of the galaxy (i.e. when 
the first star shone) and varies in the range 2 — 12 Gyrs. 
We assumed a constant small amount of extinction for all 
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populations. The relation|Sldescribes the evolution of SFR 
in an open chemical evolution model obeying the Schmidt 
law. The power-law term represents the early growth of 
SFR. 

The second method, meant to reproduce case (a), 
consisted in the same free-form approach already de- 
scribed for IR starbursts, but constrained in order to avoid 
high extinction values and high FIR emission. The lat- 
ter assumes that we have imaging data of enough quality 
to ensure that the object under scrutiny is a bona-Ede 
spheroidal galaxy not affected by extinction features as 
revealed by dust lanes, absorption bands, etc. For the typ- 
ical E/SOs in our sample, best-fits are obtained with the 
combination of just two SSPs, one of 1.2 10^" yrs and the 
other of 10^ yrs, plus a younger SSP needed for the blue 
ellipticals. 

Particularly to fit the case (b) sub-population, we have 
considered a third choice: a combination of two old single 
stellar populations, both of them 1.2 10^° yrs old, having 
different values of metallicity (solar Z=0.02 and sub-solar 
Z=0.008). 

As a final, check we have also considered the most 
general case of a completely free-form approach, as that 
mentioned above to fit the dusty star-forming class, but 
with an upper limit imposed on the mid-IR flux (as might 
be provided by forthcoming SIRTF multiband surveys) to 
constrain to amount of dust re-radiation and optical ex- 
tinction. 

3.4. Fitting the Ly-break galaxy SEDs 

The spectral energy distributions of the Ly-break galax- 
ies have been fitted assuming a discrete number of stellar 
populations as for the IR starbursts, but since at redshift 
z ~ 2.5 the age of the universe for the assumed cosmology 
is roughly 2.5 Gyrs, only SSPs up to 1 or 2 10^ years — 
depending of z — have been considered. 

As for the starburst case, only stellar populations con- 
tributing to the visible and IR light by more than 1% were 
involved in the fit. The typical number of SSPs used is 2—3 
for Lyman-break galaxies (see Table OJ. 

3.5. Parameter space exploration with the Adaptive 
Simulated Annealing (ASA ) 

For each galaxy the code randomly generates a wide set 
of models and compares them to the observed fluxes, by 
minimizing the Merit Function, which tests the accuracy 
of the given fit: 

^ / Qobs Qmod \ 

where 5'°''* and 5'™°'' are the observed and synthetic fluxes 
in the z-th photometric band, af"^ is the uncertainty on 
the observed flux and N the number of constraints. 

In the most general case there will not exist a unique 
minimum of the Merit Function, but rather a large 
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number of relative minima. The exploration of the A'^- 
dimensional parameter space has been performed adopt- 
ing the Adaptive Simulated Annealing algorithm (ASA, 
Ingber 1989, 2001), which is very well suited to search for 
a solution as close as possible to the absolute minimum, 
avoiding local minima. 

The heart of the Simulated Annealing method (Ingber 
et al., 1989, 2001) is the thermodynamic analogy with 
cooling metals: at high temperatures each molecule moves 
with respect to the others, but when T slowly decreases 
its mobility is gradually lost and matter reaches a pure 
crystalline configuration. If the cooling process is slow 
enough, metals can spontaneously reach this minimum en- 
ergy state. 

The function to be optimized can be considered similar 
to the total energy of the system. The Boltzmann proba- 
bility distribution 

Pb - e-* (7) 

shows that even at thermal equilibrium - with tempera- 
ture T - a non-zero probability to find the system in an 
energy state higher than KT does exist. Consequently it 
is possible for the system to move from the local energetic 
minimum to an excited state and subsequently evolve to 
an energy smaller than at start. 

Following this analogy, the optimization process starts 
from a given state in the A^-dimentional parameter space 
and then evolves from Ei to E2\ the probability of this 
transition is P12 — e kt . U Ei > E2, then P > 1: the 
system always moves towards lower energy states, even if 
a non-zero probability to increase its energy does exist. 
On the other hand, small values of the temperature and 
large energetic increments imply low probabilities for the 
system to evolve uphill. 

In order to converge, at each i-th iteration, the tem- 
perature of the system is adapted according to the relation 
r(i + l) = Rrp. T{i). The adopted value for Rt is 0.85, as 
suggested by Corana et al. (1987). As T declines, moves 
to excited energy states are suppressed and the algorithm 
focuses on the most promising area for optimization. On 
the other hand, when a new deep minimum is found, these 
ranges are increased, in order to avoid the system to lag 
on it. 

Minimization ends when the value decreases by 
less than the Tolerance set by the user. Further informa- 
tions on the ASA algorithm can be found at the URL 
http : //www . ingber . com/ 



4. Results 

For each of the sample sources, our simulations generate 
a very large number of models, each one based on a dif- 
ferent set of physical parameters. For a typical number of 
computed solutions 100,000), the exploration takes 
between -20 and 45 minutes on a 1.7GHz CPU PC - de- 
pending on the total number of SSPs, the adopted SFH 
and convergence criteria. Acceptable solutions (with re- 
duced Xo ~ 1) found for most of the sample galaxies. 
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Fig. 5. SEDs fits of some Lyman-break galaxies. Left: best fit solutions for sources Ly2, Ly3 and Ly5, as obtained 
involving only young SSPs. Right: consequences of introducing the > 10^ yrs old populations on the fit of source 
Lyl. The relative contribution of old stars to the total assembly of stellar mass increases from top to bottom: the fit 
accuracy does not change, but in the restframe 1 — 2 ^m domain the model becomes significatively steeper (see also 
right panel of figure [7| and text for details). 



In the case of unacceptable fits, it may be that the uncer- 
tainties in the measured fluxes were underestimated, or 
that the distribution of the errors on parameters is not 
gaussian (Papovich et al., 2001). Finally, the SSP spectra 
may have intrinsic modellistic uncertainties. 

4.1. Mid-1 R selected starburst galaxies 

Figure n illustrates the comparison of the synthetic SEDs 
to the observed datapoints and plots of the main physical 
parameters, for a representative galaxy. The left panel in 
the flgure shows the best-fits to the SEDs. The solid line is 
the sum of contributions from all SSPs; the dashed line is 
the sum of all the young SSPs (age < 10^ yrs), while the 
dotted line represents the contribution of older popula- 
tions (age > 10^ yrs); the dashed thick line longwards of 5 
/im is the redshifted M82 template; the square black dots 
are the observed photometric datapoints, the red triangle 
at A ~ 15 fim is the fit to the ISO LW3 datapoint. The 
identification number of the model (for the given source) 
and the value based on the optical/NIR fit only (LW3 
datapoint excluded) are also reported in the plots. 

The right panels of Fig. ^ show — for the various 
SSPs making up the best-fit model, whose age intervals 
are indicated in the x-axis — trends of the SFR (top his- 



tograms, y-scale on the left end), E(B-V) (top panel, filled 
squares, y-scale on the right end), the ratio of the FIR 
flux contributed by the SSP to the total flux in the V 
band [Ffir^ssp / Fv,tot, middle panel), and the cumula- 
tive distribution of the stellar mass as a function of SSP's 
age (bottom panel). The graph in the middle panel, in par- 
ticular, details the contribution of each SSPs to the FIR 
bolometric emission of the galaxy, taking into account the 
effect of dust extinction in the optical light. The bottom 
plot shows that, typically, the youngest stars contribute to 
the bulk of the FIR bolometric emission of these sources, 
but provide almost no contribution to the optically lu- 
minous mass in the galaxy. Results for additional dusty 
starbursts in Hubble Deep Field South may be found in 
Franceschini et al. (2003). 

4.2. Degeneracies in the mass estimates for IR-selected 
dusty starbursts 

The left panel of Figure El shows plots of contours as 
a function of stellar mass M and FIR luminosity LpiR, 
based on fits to the optical/NIR data only (UBVRIJHK 
photometry). The contours correspond to 1, 2, 3cr un- 
certainties. Note that even the darker la contours are ex- 
tended, showing substantial degeneracies in the solutions. 
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The horizontal lines in Fig. |21 (left panel) mark the 
measured FIR bolometric (8—1000 fim) luminosity and its 
uncertainty based on the fits to the observed LW3 flux. In 
only a fraction of the high-z galaxies the FIR constraint is 
effective in reducing the uncertainties in the stellar mass: 
this happens in particular for the least active galaxies 
(those with the lowest mid-IR excess and SFR compared 
with the optical SED). Source S27 (z = 0.58) is one such 
case: imposing the LW3 flux constraint, the stellar mass 
range is reduced by more than a factor of two. 

On the other hand S40, an ultraluminous galaxy at 
z=1.27, is a typical case for which the FIR information 
is barely useful to constrain the mass. This fact is due 
to the presence of highly obscured young stellar popula- 
tions in the most active starbursts, which give a significant 
contribution to the NIR and FIR luminosity, but are not 
visible in the optical/NIR, because of their low M/L ratio, 
and their contribution to the stellar mass budget is hard 
to constrain. For such complex situations, further mid-IR 
data from SIRTF will be critical, as discussed in Sect. |S1 
below. 

As mentioned in Sect. 3.2, we have reproduced the con- 
tribution by old stellar populations with spectra obtained 
assuming constant SF between 2 and 12 Gyrs, since the 
single SSP spectra in this age range are all very similar 
and their M/L ratio varies smoothly. This allows us to re- 
duce the number of free parameters, when compared to a 
more detailed SFH. To test this procedure, we have ran 
different simulations on several sources, detailing the SFH 
between 2 and 12 Gyrs into 5 different SSPs (2, 4, 7, 10 
and 12 Gyrs old). Other simulations obtained by simply 
substituting the long-SF model with single 2 and 12 Gyrs 
SSPs have been run as well, in order to test the more 
extreme cases of masses assembled with the lowest and 
the highest M/L values in the 2-12 Gyrs range. The most 
massive and the lightest models obtained in this way may 
differ by a factor of ^ 2 in mass. However, this analysis 
showed that the uncertainty AM in stellar mass - due to 
different adopted models of the old population - is in fact 
completely contained within those of the global model fit. 
Indeed the AM due to the presence of young stellar pop- 
ulations dominates the global mass uncertainty. 

Although we have based our analysis on solar- 
metallicity SSPs, we have performed a systematic anal- 
ysis of how variations in metallicity may affect our re- 
sults, taking into account that e.g. changing it from so- 
lar to one third solar would be equivalent to change the 
SSP age from 10^ to 12 x 10^ yrs (the age-metallicity de- 
generacy). In principle one would expect a corresponding 
significant difference in the M/L ratio, hence on the esti- 
mated stellar mass. In fact our analysis has shown that a 
kind of compensations holds, such that changing the stel- 
lar metallicity requires a substantial modification of the 
SF history to keep an acceptable spectral fit. The result 
is that the mass estimate does not significantly depend 
on the assumed metal content, and the same applies for 
the mass uncertainty range. For star-forming galaxies the 



overwhelming uncertainty factor is instead the incidence 
of the young stellar population. 

The stellar mass estimates are reported in Tabled We 
plot in the right panel of Fig.|21the values as a function 
of mass for acceptable models, each dot representing a 
different solution. Columns 8 and 9 of Table report the 
estimated uncertainties AM at 1 and 2<t levels. The typical 
uncertainty in the mass values for these IR starbursts is 
a factor of 2—3, but may be as high as 5 in a few cases. 
Figure El shows the estimated masses against redshift. 

4.3. Elliptical galaxies 

Franceschini et al. (1998) and Rodighiero et al. 
(2001) found that populations of morphologically-selected 
spheroidal galaxies in the HDFs are characterized by a 
wide variety of integrated colors, including typical ellip- 
ticals with very red SEDs (B-I>2.5 for objects at z < 1 
and V-J>2.5 for z > 1), but also many bluer objects. The 
b/r fiag in the third column of Table |21 identifies blue and 
red objects in our sample. 

As previously mentioned, the SEDs of elliptical and SO 
galaxies have been fit alternatively by assuming an ana- 
lytic star formation history, or by means of the free-form 
approach, or the sub-solar metallicity model. The SEDs of 
the red subsample are equally well-fitted by all three, and 
the corresponding estimated values of the stellar masses 
agree very well within the relatively small uncertainties. 
For these sources, the mass range reported in Table El in- 
cludes the results of the three fits, while the best fit value 
refers to the analytic SFH tecnique. 

Figure O shows best fits for the three methods for 
source E5: all three yield equally good fits, although the 
analytic SFH gives the lowest x^; the lower panels of Fig. 

report the star formation histories and the cumulative 
mass distribution in the three cases. Table 12 summarizes 
our results; the estimated stellar masses are in fair agree- 
ment with those by Rodighiero et al. (2001). 

For the blue E/SO galaxy population, the simple com- 
bination of 2 old populations cannot reproduce the ob- 
served U, B (and sometimes V) observed fluxes, while 
good solutions can be found by assuming some ongoing 
star formation (either within the analytic model, see Fig. 

01 or by adding a third young SSP in the discrete SFH). 

Finally we have investigated the implications of 
analysing the SEDs for this class of galaxies with the 
most general free-form spectro-photometric tool based on 
the multi-SSP approach and including stronger extinction 
effects. The only assumption in this experiment is that 
we have an upper limit to the mid- and far-IR flux to 
constrain possible dust re-radiation (we put a limit of 70 
/iJy at 24 /im, a value easily achievable with SIRTF ob- 
servations). For red ellipticals this increased freedom in 
the spectral analysis has no effect in modifying the stellar 
mass estimate, while for blue ellipticals the uncertainty 
range increases by typically 20%, with respect to the case 
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Fig. 6. Comparison of stellar masses between mid-IR selected starburst (filled squares), K-band selected E/SO (open 
circles) and Ly-break galaxies (triangles), as a function of redshift. With the exception of very massive MIR sources 
belonging to the HDFS well-known group at z ~0.6, the stellar masses we find for starbursts and E/SO are comparable; 
the small sample of Ly-break galaxies seems to consist of less massive objects than those observed at lower redshift. 



when we can exclude a-priory important extinction e.g. 
based on imaging data. 

4.4. Comparison of stellar masses for moderate-z 

starbursts and E/SO and high-z Ly-break galaxies 

Figure [S] reports the best fit solutions for four of the 
Lyman-break sample galaxies based on models including 
only populations younger than 10^ yrs. The right-hand 
panels, top to bottom, of the figure illustrate the conse- 
quence of introducing an increasing fraction of 1-2 Gyr-old 
SSP in the spectrum of source Lyl. The overall fit to the 
observed optical-NIR data and values do not change, 
but the synthetic SED longwards of the K-band is signif- 
icantly modified. As a consequence, without constraints 
on the rest-frame near-IR spectrum where the bulk of the 
emission by low-mass stars is expected, the uncertainties 
in the photometric mass estimates are very large, typically 
1 order of magnitude (see numerical details in Table 

Figure plots the estimated baryonic masses for our 
three galaxy samples as a function of redshift. On one side, 
this provides evidence that the IR-selected starbursts and 
the K-band selected ellipticals and SO galaxies seem to be 
comparably massive systems. Note however that the small 
sampled volume may imply important cosmic variance ef- 
fects (e.g. some structures are evident in the HDF South 
samples around z=0.5-0.6, including galaxies of huge stel- 
lar masses). 



Fig.|H|shows that, in spite of the large mass uncertain- 
ties, our small illustrative sample of Lyman-break galax- 
ies seem to correspond to a galaxy population signifi- 
cantly less massive {M < a few IO^'^Mq) then those ob- 
served at lower redshifts (for which typically M > several 
IO^^Mq, except for z<0.7 where low-mass galaxies be- 
come detectable) over comparable sky areas. Similar re- 
sults were found by Papovich et al. (2001) on the sample 
of HDFN faint Lyman-break galaxies, with a median stel- 
lar mass of M ^ 6 10^ Mq and by Shapley et al. (2001) on 
brighter Ly-break galaxies (with a slightly higher median 
mass). Although not conclusive, due to the lack of statis- 
tics and completeness, these results indicate substantial 
stellar build-up to happen at z ~ 1 to 2 in the field pop- 
ulation. This seems consistent with the steady increase of 
galaxy's M/L from z=3 to z=0 and the comoving stellar 
mass density increasing by a factor of ^ 10 between z=2 
and z=l claimed by Dickinson et al. (2003). 

5. Perspectives for the SI RTF mission 

The newly operative NASA's Space InfraRed Telescope 
Facility (SIRTF) includes a sensitive camera (IRAC) to 
observe high-z galaxies in complementary long- wavelength 
channels between 3 and 10 /im, hence providing us with 
a unique opportunity for the analysis of stellar masses at 
intermediate and high (0.5—3) redshifts, by sampling the 
critical restframe near-IR domain. 
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Fig. 7. Simulations of SIRTF/IRAC observations to constrain the estimate of stellar masses in starbursts and Lyman- 
break galaxies. Left panel: three different solutions for IS0/LW3 source S40: all three fits reproduce comparatively 
well both the optical data and the LW3 observed flux (even if not displayed for plotting reasons). SIRTF synthetic 
data (open squares) have been obtained simulating an IRAC observation of the best fit model. Right panel: the three 
different fits to optical- NIR data of source Lyl already discussed in Fig. IRAC synthetic fluxes have been computed 
for the best fit solution, which is the youngest and least massive among the three considered. In both panels the three 
models differ in their star formation history: the reddest are dominated by old populations more than the others. Small 
boxes show the transmissivity curves of the four IRAC filters; the numbers of the models refer to those appearing in 
Figs. H] and 



We have simulated SIRTF/IRAC observations based 
on the best fit models of some sources in the selected 
starburst and Lyman- break galaxy samples, to attempt 
to quantify how such observations will improve the mass 
estimates. 

5.1. Constraining the mass of active starbursts at 
z = 0.5- 1.5 

In the case of the IR-selected starbursts, we have found 
uncertainties in the stellar mass estimate based on 
the optical/NIR+LW3 SED of typically factors of 3-4. 
SIRTF/IRAC, observing at 3.6, 4.5, 5.8 and 8.0 ^m, will 
significantly improve upon this. Figure \7\ illustrates this 
with an example. The open squares are the expected fluxes 
for galaxy S40 in the SIRTF/IRAC bands, assuming our 
best-fit spectral solution. Some other solutions, character- 
ized by different SFHs, are also plotted. The case of S40 
is representative of the high redshift population, while for 
intermediate z objects (e.g. S53 or S55, z ^ 0.5 — 0.7) the 
improvement may be less significant. 

We have simulated a SIRTF/IRAC observation of 
twelve randomly-selected sources in the LW3 dusty star- 
burst sample. We have convolved the best-fit optical/NIR 
SED with the 4 IRAC band transmission curves. Then we 
have run our code of spectral synthesis on the 12 objects, 
assuming IRAC simulated fluxes with a 10% photometric 



errors. As a result, the estimated stellar-mass range is ex- 
pected to reduce to a typical factor of ~ 2, as reported in 
the last column of Table ^ 



5.2. Detecting NIR restframe emission of z = 2 — 3 
Lyman-break galaxies 

Lying at redshifts z > 2, the Lyman-break galaxies in our 
sample will benefit by IRAC observations much more than 
intermediate and low-z galaxies. In this redshift range, 
IRAC 3 — 8 /im data will sample the restframe near-IR 
emission of galaxies, hence observing the light emitted by 
low-mass stars dominating the galaxy stellar content. 

We simulated IRAC observations based on our best 
fit solutions. Then we have applied our spectral synthe- 
sis code using 11 datapoints, namely the already available 
UBVIJHK photometric data and the synthetic IRAC 3.6, 
4.5, 5.8 and 8.0 fim fluxes. Figure [3 (right-hand panel) 
shows how IRAC observations will constrain the three so- 
lutions shown in the right panel of figure |S1 

The results of these simulations are summarized in the 
last column of Table |21 the mass uncertainty very signifi- 
cantly reduces to a factor of ~ 2 — 3. 



12 



Berta S., et al.: Photometric mass estimates in high-z galaxies 



6. Summary 

We have presented a new spectrophotometric synthesis 
code for the integrated Hght of distant and high-redshift 
galaxies. This tool is aimed, in perspective, at investigat- 
ing various fundamental physical parameters, like the age 
of the dominant stellar populations, extinction, SF history, 
ongoing star-formation rate, etc. We have concentrated in 
this paper to study in particular the ability of the code 
to constrain the stellar mass from SEDs fitting. The tool 
has been tested on a sample of luminous infrared star- 
bursts detected by ISO in the Hubble Deep Field South 
at 15 /xm (Franceschini et al., 2003), on a HDFS K-band 
selected sample of ellipticals and SOs (Rodighiero et al. 
2001), and on a set of the Lyman-break galaxies in the 
HDF South from Papovich et al. (2001). 

We have fitted the optical-NIR SEDs (from U to K 
bands) of the mid-IR starbursts by assuming a discrete 
star formation history between 1.2 10^" years ago and to- 
day, by means of the combination of up to 10 single stellar 
populations. Each SSP was weighted by different SFR's 
and absorbed by different amounts of dust. For a best ex- 
ploitation of the available constraints we also included in 
the analysis the 15-/im ISO flux, to assess the contribu- 
tion of extinguished young stars to the emitted spectrum. 
Based on these data, our estimated mass uncertainties 
range from a factor ~2 to occasionally a factor of 5. In such 
dust-obscured starburst galaxies, substantial uncertainty 
in the mass estimate comes from the still possible exis- 
tence of strongly extinguished young stars with low M/L 
ratios, contributing to the far-IR flux but undetectable in 
the rest-frame optical. 

We have also analysed a sample of morphologically- 
selected E/SO galaxies, at z<1.5. Their SEDs have been 
reproduced as the combination of stellar populations with 
solar metallicities, or by combining SSPs of different 
metallicities (Z=0.02 and Z=0.008). Additional solutions 
have been sought with a parameterized continuous se- 
quence of stellar ages. The SEDs of spheroidal galaxies 
with red colors are well fitted by combining intermedi- 
ate and old populations only, while for ellipticals with 
relatively blue colors (the majority in our complete sam- 
ple) acceptable solutions are found only by assuming some 
some recent or even ongoing star formation. The three dif- 
ferent methods provide results well consistent with each 
other. 

For Lyman-break galaxies the available rest-frame 
optical-UV data are much less constraining. The conse- 
quences of assuming relatively old 10^ yrs) stellar pop- 
ulations in these systems have been discussed and it has 
been found that the lack of long-wavelength data imply 
uncertainties in the mass estimate of up to a factor of 

10. In spite of these large uncertainties, there are indi- 
cations that typical Lyman-break galaxies correspond to 
a galaxy population significantly less massive then those 
observed at lower redshifts, possibly indicating substan- 
tial stellar build-up to ocurring at z ^ 1 to 2 in the field 
galaxy population. 



The situation is expected to significantly improve, 
particularly for higher-z galaxies, with the forthcoming 
SIRTF/IRAC observations. Our simulations show that 
such an extension of the wavelengths dynamic range will 
reduce the mass uncertainties to factors of 2— 3 for var- 
ious classes of galaxies up to z~3. This promises to be 
good enough for statistically reliable determinations of the 
galaxy evolutionary mass functions. 
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^ 09 


0.016-0.044 


0.013-0.067 






S72 


0.55 


421227 


3 


1.40 


12.39 


1.14-1.30 


1.12-2.01 






S73 


0.17 


441608 


4 


1.50 


10.00 


0.95-1.01 


0.84-1.94 


0.96- 


-1.82 


S75 


(0.45) 


222010 


4 


1.30 


4.19 


1.16-1.36 


0.83-1.40 






S77 


(0.40) 


210033 


4 


0.80 


4.73 


0.91-1.10 


0.47-1.16 


0.45- 


-1.14 


S79 


0.74 


499221 


3 


0.60 


1.60 


0.65-1.18 


0.30-1.21 






S82 


0.69 


438044 


4 


0.50 


9.61 


0.49-0.62 


0.37-0.83 


0.51- 


-0.69 


S85 


(0.40) 


414053 


3 


5.70 


25.77 


5.51-5.79 


4.78-6.19 







Table 1. Results of the stellar-mass estimate based of SEDs fitting. For each IS0/LW3 HDFS source we report redshift 

(Franccschini ct al. 2003; photometric if within parenthesis, spectroscopic otherwise), total number of models produced 
by the spectrophotometric synthesis code (without the SIRTF canstraints), effective number of SSPs contributing to 
best fit (for the optical -I-LW3 data, but without the SIRTF constraints), optical and mass of the best fit solution. 
The last three columns report three different estimates of the mass uncertainties due to degeneracies: AM at 68.3% 
and 95.5% confidence levels (corresponding to la and 2a) and 2a mass range obtained simulating SIRTF/IRAC 
observations (only for some randomly selected sources). 
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Obj 


z 


Color 


Ntot 


tF 


E(B-V) 


ni 


Tl 




Mass(bf) 


AM 


# 




Flag 


mdls. 


b.f. [Gyrs] 


b.f. 


b.f. 


b.f. 


b.f 


[10" M©] 


[10" M©] 


El 


0.512 


b 


18001 


11.91 


0.05 


2.69 


0.16 


8.67 


0.18 


0.08-0.18 


E2 


(1.3) 


b 


16848 


11.66 


0.02 


3.01 


0.16 


11.95 


0.83 


0.47-1.07 


E3 


(1.3) 


b 


23427 


3.44 


0.05 


4.12 


0.08 


7.95 


0.85 


0.70-1.82 


E4 


0.414 


b 


11521 


2.68 


0.00 


0.002 


0.52 


4.29 


0.16 


0.10-0.23 


E5 


0.673 


r 


129606 


3.19 


0.05 


0.001 


0.15 


1.89 


0.73 


0.66-1.12 


E6 


(0.5) 


r 


111765 


9.20 


0.06 


0.002 


0.15 


6.31 


0.90 


0.52-1.16 


E7 


0.517 


b 


11287 


4.63 


0.10 


5.17 


0.11 


6.71 


0.10 


0.08-0.14 


E8 


0.579 


r 


131529 


6.68 


0.05 


2.85 


0.09 


5.31 


0.52 


0.42-0.83 


E9 


0.515 


b 


15841 


2.00 


0.00 


5.25 


0.10 


6.68 


0.05 


0.045-0.059 


Ell 


0.364 


r 


121526 


7.17 


0.04 


2.35 


0.06 


9.92 


0.10 


0.07-0.16 


E14 


(1.15) 


b 


11053 


11.44 


0.00 


1.97 


0.25 


5.52 


0.75 


0.43-0.81 


E16 


(1.0) 


b 


11283 


11.46 


0.05 


0.04 


0.26 


8.36 


3.42 


2.29-3.99 


E17 


(1.25) 


r 


136281 


11.71 


0.04 


0.03 


0.19 


7.31 


3.89 


2.82-4.08 


E18 


(0.9) 


r 


112240 


11.46 


0.04 


0.02 


0.19 


4.58 


0.73 


0.38-0.89 


E19 


(0.15) 


r 


121521 


11.96 


0.05 


3.98 


0.01 


9.81 


0.08 


0.06-0.09 


E20 


(0.8) 


b 


16087 


2.00 


0.00 


0.001 


0.34 


9.76 


0.12 


0.08-0.18 


E21 


0.577 


r 


127454 


11.28 


0.05 


0.04 


0.29 


10.00 


0.98 


0.61-1.16 


E24 


(0.4) 


r 


134001 


11.28 


0.05 


3.35 


0.08 


13.54 


0.55 


0.33-0.72 


E25 


0.511 


b 


13272 


11.97 


0.04 


1.49 


0.21 


7.02 


0.72 


0.44-0.93 


E26 


(1.5) 


b 


11729 


5.52 


0.00 


3.69 


0.09 


22.34 


2.80 


1.89-4.12 


E27 


(0.85) 


b 


11258 


11.72 


0.02 


3.42 


0.11 


16.02 


0.92 


0.52-1.09 


E28 


0.564 


b 


11281 


11.91 


0.02 


0.64 


1.70 


18.45 


0.13 


0.086-0.14 


E29 


(0.5) 


r 


112123 


14.35 


0.05 


2.76 


0.01 


4.79 


0.18 


0.078-0.32 



Table 2. Results of the stellar-mass estimate for the ellipticals sample, based of SEDs fitting with solar metallicity 
SSPs. For each elliptical source in the HDFS (Rodighiero et al., 2001), spectroscopic (Vanzella et al. 2003, Sawicky & 
Mallen-Ornelas 2003) or photometric (if within parenthesis, Rodighiero) redshifts, total number of models produced by 
the spectrophotometric synthesis code, best fit parameters and for the analytic SFH method (see text for details), 
mass of the best fit solution and 2a (95.5% confidence) mass range due to degeneracies are reported. Third column 
reports a flag identifying red and blue sources. The SEDs of the red sources have been fitted with three different 
methods: a) analytic SFH, h) two single stellar populations (10^ and 1.2 10^*^ yrs old) with solar metallicity, c) two old 
populations (1.2 10^° yrs) with different metallicities (Z=0.02, 0.008) — see also text and figure|3| The mass ranges 
reported for these objects include the results of all three fitting tecniques; on blue sources only the analytic SFH leads 
to good results. 



Obj 


ID Pap. 


ID Will. 


z 


Ntot 


Nssp 


'2 

Xopt 


Mass(bf) 


AM 


AMsiRTF 


# 


(1) 


(2) 


(3) 


mdls. 


b.f. 


b.f. 


[10" Mq] 


[10" Mq] 


[1O"M0] 


Lyl 


110 


2-449.0 


2.005 


230403 


3 


5.13 


0.21 


0.16-2.47 


0.27-0.55 


Ly2 


503 


2-903.0 


2.233 


195602 


3 


4.24 


0.13 


0.04-0.24 


0.09-0.19 


Ly3 


67 


2-82.1 


2.267 


195587 


3 


1.47 


0.04 


0.03-0.84 


0.04-0.09 


Ly4 


804 


4-639.0 


2.591 


212778 


2 


10.99 


0.09 


0.03-0.28 


0.03-0.09 


Ly5 


1352 


4-497.0 


2.800 


199203 


3 


5.87 


0.02 


0.02-0.29 


0.02-0.06 


Ly6 


1541 


4-363.0 


2.980 


217529 


3 


9.02 


0.18 


0.05-0.35 


0.06-0.20 


Ly7 


661 


2-643.0 


2.991 


264435 


2 


7.54 


0.04 


0.02-0.13 


0.03-0.06 


Ly8 


273 


2-76.11 


3.160 


278399 


2 


2.29 


0.56 


0.27-1.35 


0.38-0.81 



(1) : ID number in Papovich et al. (2001); 

(2) : ID number in Williams et al. (1996); 

(3) : Spectroscopic redshifts by Steidel et al (1996) and Lowental et al. (1997). 



Table 3. Results of the spectral synthesis analysis of the lyman-break galaxies. Objects are identified by numbers 
adopted in this work, by Papovich and by Williams (first three columns). For each source the table reports redshift, 
total number of models analyzed, effective number of SSPs involved in the best fit, stellar mass and of best fit and 
2a mass range due to degeneracies. Last column summarizes the results of our simulated SIRTF/IRAC observations 
of ly-break sources: IRAC will sample the restframe near-IR light at z 2 — 3, therefore significantly reducing the 
uncertainty on the stellar mass estimate (see text for details). 
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